d-wave superconductivity and Pomeranchuk instability 
in the two-dimensional Hubbard model 
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We present a systematic stability analysis for the two-dimensional Hubbard model, which is 
based on a new renormalization group method for interacting Fermi systems. The flow of effective 
interactions and susceptibilities conflrms the expected existence of a d-wave pairing instability driven 
by antiferromagnetic spin fluctuations. More unexpectedly, we find that strong forward scattering 
interactions develop which may lead to a Pomeranchuk instability breaking the tetragonal symmetry 
of the Fermi surface. 
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The two-dimensional Hubbard model |i| has attracted 
much interest as a promising prototype model for the 
electronic degrees of freedom in the copper-oxide planes 
of high-temperature superconductors, since it has an 
antiferromagnetically ordered ground state at half-filling 
and is expected to become a d-wave superconductor for 
slightly smaller electron concentrations 

Although the Coulomb interaction in the cuprate su- 
perconductors is rather strong, the tendency towards an- 
tiferromagnetism and d-wave pairing is captured already 
by the 2D Hubbard model at weak coupling. Conven- 
tional perturbation theory breaks down for densities close 
to half-filling, where competing infrared divergences ap- 
pear as a consequence of Fermi surface nesting and van 
Hove singularities [^-D . A controlled and unbiased treat- 
ment of these divergencies cannot be achieved by stan- 
dard resummations of Feynman diagrams, but requires 
a renormalization group (RG) analysis which takes into 
account the particle-particle and particle-hole channels 
on an equal footing. 

Early RG studies of the two-dimensional Hubbard 
model started with simple but ingenious scaling ap- 
proaches, very shortly after the discovery of high-Tj, su- 
perconductivity ^-|^. These studies focussed on dom- 
inant scattering processes between van Hove points in 
k-space, for which a small number of running couplings 
could be defined and computed on 1-loop level. Spin- 
density and superconducting instabilities where identi- 
fied from divergencies of the corresponding correlation 
functions. 

A major complication in two-dimensional systems 
compared to one dimension is that the effective inter- 
actions cannot be parameterized accurately by a small 
number of running couplings, even if irrelevant momen- 
tum and energy dependences are neglected, since the tan- 
gential momentum dependence of effective interactions 
along the Fermi surface is strong and important in the 
low-energy limit. This has been demonstrated in par- 
ticular in a 1-loop RG study for a model system with 
two parallel flat Fermi surface pieces ||^. Zanchi and 
Schulz H have recently shown how modern functional 



renormalization group methods can be used to treat the 
full tangential momentum dependence of effective inter- 
actions for arbitrary curved Fermi surfaces. Most re- 
cently, Salmhofer ||^ has derived an improved version of 
this field theoretic approach. The resulting flow equa- 
tions are particularly suitable for a concrete numerical 
evaluation. To compute physical instabilities, we have 
derived the corresponding flow equations for susceptibil- 
ities §. 

In this letter we present results for the flow of suscep- 
tibilities as obtained by applying Salmhofer 's renormal- 
ization group method to the two-dimensional Hubbard 
model with nearest and next-nearest neighbor hopping 
on a square lattice. The expected existence of a d-wave 
pairing instability driven by antiferromagnetic spin fluc- 
tuations is thereby conflrmed beyond doubt. More un- 
expectedly, we find that strong forward scattering inter- 
actions develop which may lead to a Pomeranchuk |Tc| ] 
instability breaking the tetragonal symmetry of the Fermi 
surface. 

The one-band Hubbard model 101 
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describes tight-binding electrons with a local repulsion 
[/ > 0. Here c\^ and cia are creation and annihilation 
operators for fermions with spin projection cr S {t, i} on 



a lattice site i, and 
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A hopping amplitude 



-t between nearest neighbors and an amplitude —t' be- 
tween next-nearest neighbors on a square lattice leads to 
the dispersion relation 

Ck = — 2i(cos kx + cos ky) — 4t' cos k^ cos ky (2) 

for single-particle states. This dispersion relation has 
saddle points at k = (0,7r) and (tt, 0), which generate 
logarithmic van Hove singularities in the non-interacting 
density of states at the energy e^H = 4t'. For t' = 0, 
Ek has the nesting property ek+q — — Ck for Q = (tTjTt), 
which leads to an antiferromagnetic instability for arbi- 
trarily small [/ > at half-filling . 
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The RG equations are obtained as follows (for details, 
see Salmhofer [|| and Ref. ||]). The infrared singularities 
are regularized by introducing an infrared cutoff A > 
into the bare propagator such that contributions from 
momenta with jet — /i| < A arc suppressed. All Green 
functions of the interacting system will then flow as a 
function of A, and the true theory is recovered in the 
limit A ^ 0. Salmhofer [|j has recently pointed out that 
(amputated) Green functions obtained by expanding the 
effective action of the theory in powers of normal ordered 
monomials of fermion fields obey differential flow equa- 
tions with a structure that is particularly convenient for 
a power counting analysis to arbitrary loop order. With 
the bare interaction as initial condition at the highest 
scale Ao — maxjet — these flow equations determine 
the exact flow of the effective interactions as A sweeps 
over the entire Brillouin zone down to the Fermi surface. 
The effective low-energy theory can thus be computed 
directly from the microscopic model without introducing 
any ad hoc parameters. 

For a weak coupling stability analysis it is sufficient 
to truncate the exact hierarchy of flow equations at 1- 
loop level. The effective 2-particle interaction then re- 
duces to the one- particle irreducible 2-particle vertex F^, 
and its flow is determined exclusively by F"^ itself (no 
higher many-particle interactions enter) . Flow equations 
for susceptibilities are obtained by considering the exact 
RG equations in the presence of suitable external fields, 
which leads to an additional 1-particle term in the bare 
interaction, and expanding everything in powers of the 
external fields to sufficiently high order |^ . 

One cannot solve the flow equations with the full en- 
ergy and momentum dependence of the vertex function, 
since F^ has three independent energy and momentum 
variables. The problem can however be much simpli- 
fied by ignoring dependences which are irrelevant in the 
low energy limit, namely the energy dependence and the 
momentum dependence normal to the Fermi surface (for 
details, see Ref. |^). This approximation is exact for the 
bare Hubbard vertex, and asymptotically exact in the 
low-energy regime. The remaining tangential momen- 
tum dependence is discretized for a numerical evaluation. 
Most of our results where obtained for a discretization 
with 16 points on the Fermi surface (yielding 880 "run- 
ning couplings" ) , and we have checked that increasing the 
number of points does not change our results too much. 

We have computed the flow of the vertex function for 
many different model parameters t' and U (t just fixes the 
absolute energy scale) and densities close to half-filling. 
In all cases the vertex function develops a strong momen- 
tum dependence for small A with divergencies for several 
momenta at some critical scale Ac > 0, which vanishes 
exponentially for U 0. To see which physical insta- 
bility is associated with the diverging vertex function we 
have computed commensurate and incommensurate spin 
susceptibilities xs(q) with q = (tt, tt), q = (tt — 6, tt) and 
q = (1 — i5)(7r,7r), where S is afunction of density the 
commensurate charge susceptibility xc(7i'j "■), and singlet 
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FIG. 1. The flow of the ratio 
non-interacting susceptibilities for t' = 
M = -0.055t. 



of interacting and 
-O.Olt, [/ = t and 



pair susceptibilities with form factors 



sin kr sin fc,, 



i (s-wave) 
I -^{coa kx + cos ky) (extended s-wave) 
"^^^^ I -^(cos fca; — cos fcjy) (d-wave (ij,2_j,2) 

(d-wave dxy). 
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Some of these susceptibilities diverge together with the 
vertex function at the scale A^,. Depending on the choice 
of U, t' and fj,, the strongest divergence is found for the 
commensurate or incommensurate spin susceptibility or 
for the pair susceptibility with dx^^y2 symmetry. In Fig. 
1 we show a typical result for the fiow of susceptibili- 
ties as a function of A. Note the threshold at A w 0.03f 
below which the amplitudes for various scattering pro- 
cesses, especially umklapp scattering, renormalize only 
very slowly. The fiow of the antiferromagnetic spin sus- 
ceptibility is cut off at the same scale. The pairing sus- 
ceptibility with c?2,2_y2-symmetry is obviously dominant 
here (note the logarithmic scale). Following the fiow of 
the susceptibilities one can see that the (ia;2_y 2 -pairing 
correlations develop in the presence of pronounced but 
short-range antiferromagnetic spin- correlations, in agree- 
ment with earlier ideas on c?-wave superconductivity in 
the Hubbard model 

In Fig. 2 we show the (/i, U) phase diagram for t' ~ 
— G.Olt obtained by identifying the dominant instability 
from the fiow for many different values of fi and U . For 
fi = At' the Fermi surface touches the saddle points (0, tt) 
and (tt, 0), while fi = 4t'-|-0.01t corresponds to half-filling. 
Note that for [/ ^ the pairing instability always dom- 
inates, because the BCS channel dominates the fiow in 
the limit A — > 0. A spin density wave is the leading in- 
stability for U ^ only in the special case with perfect 
nesting, t' — and fi = (see the (/i, U) phase diagram 
computed from the 1-loop flow for t' = in Ref. Q). 
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FIG. 2. The (/x, [/) phase diagram for t' = -O.Olt near 
half-fiUing (marked by the dashed vertical line); the symbols 
represent the parameter values for which the flow has been 
computed and whether the dominant instability is magnetic 
(squares) or superconducting (circles); the solid line sepa- 
rates the spin-density wave regime from the superconducting 
regime. 
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FIG. 3. The critical energy scale Ac as a function of the 
chemical potential /i for U = 1.5t and t' = — O.Olt. The dif- 
ferent symbols indicate whether the leading instability is a 
spin-density wave or d-wave pairing instability. 

How the critical energy scale Ac varies as a function 
of the chemical potential (i.e. as a function of density) 
is shown in Fig. 3 for an interaction strength U ~ l.bt. 
Obviously Ac is maximal for a chemical potential at the 
van Hove energy. Note that Ac must not be interpreted 
as a transition temperature for spin density wave forma- 
tion or superconductivity, but rather as an energy scale 
where bound particle-hole or particle-particle pairs are 
formed. 

Since some of the forward scattering interactions grow 
strong for small A, while the Fermi velocity is very small 
near the saddle points, the Fermi surface may be signifi- 
cantly deformed by interactions, especially for /i « e^n- 





-n 

(a) -71 

FIG. 4. Schematic 
breaking the square symmetry; the deformed surface may be 
closed (a) or open (b). 
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plot of Fermi surface deformations 



Previous investigations of Fermi surface deformations 
within standard perturbation theory have yielded only 
very small shifts even for sizable interaction strengths 
|l2| , but in these studies the possibility of a spontaneous 
breaking of the point group symmetry of the square lat- 
tice has not been taken into account. 

To analyze systematically the stability of the Fermi 
surface shape, we define a susceptibility 



(4) 



which measures the size of Fermi surface shifts Ss^^ for 
small momentum dependent shifts of the chemical po- 
tential S/iy^i^ at points on the Fermi surface. The 
matrix Kkj^k^ defines a linear integral operator acting 
on functions of kj?. A simple consideration in the spirit 
of phenomenological Fermi liquid theory shows that the 
corresponding inverse operator is given by 



(«;"')k,k' =i^k,<^(kF-ki.) + 2/, 



kpk' 



(5) 



where v^p is the Fermi velocity and /^^j.' Landau 
function in the charge (spin-symmetric) channel. It is 
now obvious that the matrix Kki^k'^ is symmetric. The 
Fermi surface is stable, if all eigenvalues of k (or k~^) 
are positive. Note that Landau's energy functional can 
be written as a quadratic form in Ss^p, with as ker- 
nel p3[ |, and negative eigenvalues would imply that this 
energy can be lowered by a suitable deformation of the 
Fermi surface. In isotropic Fermi liquids such instabil- 
ities occur for strongly negative Landau parameters, as 



first pointed out by Pomeranchuk [ 
We have computed the renorma. 



10[. 

ization group flow of 



the eigenvalues and eigenvectors of the operator n from 
the flow of the Landau function f^^^^r , which is given di- 
rectly by the vertex function in the forward scattering 
channel |0. For various choices of the model parame- 
ters we have always found a Fermi surface instability at a 
scale Af above the scale Ac where the vertex function di- 
verges. In all cases the instability favors a deformation of 
the Fermi surface which breaks the point group symme- 
try of the square lattice, as shown schematically in Fig. 
4. The instability is mainly driven by a strong attrac- 
tive interaction between particles (or holes) on opposite 
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corners of the Fermi surface near the saddle points and 
a repulsive interaction between particles on neighboring 
corners. 

The above diagnosis of Pomeranchuk instabilities 
would be rigorous for a normal Fermi liquid with finite 
renormalized interactions in the infrared limit. In the 
present system, however, the vertex function diverges at 
a finite scale and possible Pomeranchuk type instabilities 
compete with magnetic and superconducting instabili- 
ties. Since we have no quantitative theory of the strong 
coupling physics near and below the scale A^,, we can only 
list and discuss two possible scenarios: 

i) Energy gaps due to particle-particle or particle-hole 
binding may stop the flow of forward scattering interac- 
tions before a Pomeranchuk instability sets in. 

ii) The Pomeranchuk instability is not blocked by bind- 
ing phenomena. In that case one would have a finite 
temperature phase transition with a spontaneous break- 
ing of the (discrete) tetragonal symmetry of the square 
lattice, and subsequent continuous symmetry breaking 
associated with magnetic order or superconductivity in 
the ground state. 

Which of the two scenarios is realized depends on the 
choice of the model parameters. The Pomeranchuk insta- 
bility occurs more easily if the Fermi surface is close to 
the saddle points of Ck. On the other hand, nesting raises 
the scale for particle-hole binding (leading ultimately to 
magnetic order). The best candidate is therefore the 
Hubbard model with a sizable t' (reducing nesting) and 

We emphasize that the Pomeranchuk instability does 
not cut off the singularity in the Cooper channel since 
it does not break the reflection invariance. Hence, at 
sufficiently large doping away from half-filling, d-wave 
superconductivity will set in in any case, with an or- 
der parameter that may be slightly distorted away from 
perfect d-wave symmetry. The Pomeranchuk instability 
would also not destroy the umklapp scattering route to 
an insulating spin liquid discussed recently by Furukawa 
et al. 

To our knowledge a Pomeranchuk instability has not 
yet been observed in numerical solutions of the two- 
dimensional Hubbard model. Of course this may be 
due to finite size limitations or too high temperatures 
in Monte Carlo simulations. It would thus be interest- 
ing to compute the Fermi surface susceptibility Kkj^k^ 
numerically. 

In real systems a Pomeranchuk instability as in Fig. 
4 may lead to an orthorhombic lattice distortion, as a 
consequence of the coupling of electronic and lattice de- 
grees of freedom. High temperature superconductors in- 
deed exhibit structural phase transitions between tetrag- 
onal and orthorhombic phases. It would be interesting 
to clarify whether a Pomeranchuk instability might drive 
(at least to a significant extent) the transition into the 
orthorhombic phase in these materials. 



In summary, we have shown that modern renormal- 
ization group methods can be used to establish the ex- 
pected d-wave pairing instability in the two-dimensional 
Hubbard beyond doubt. Note that for small bare inter- 
actions and in a parameter regime where only particle- 
particle pairing fluctuations grow strong, the strong cou- 
pling problem associated with the formation of a super- 
conducting state can be treated rigorously ||l6|| . Further- 
more, we have pointed out that a Pomeranchuk instabil- 
ity breaking the tetragonal symmetry of the Fermi sur- 
face is likely to occur for a suitable choice of the model 
parameters, especially for a Fermi surface close to the 
saddle points in the absence of perfect nesting. 
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